Curvature effect on nuclear "pasta": Is it helpful for gyroid appearance? 
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In supernova cores and neutron star crusts, nuclei are thought to deform to rodlike and slablike 
shapes, which are often called nuclear pasta. We study the equilibrium properties of the nuclear 
pasta by using a liquid drop model with curvature corrections. It is confirmed that the curvature 
effect acts to lower the transition densities between different shapes. We also examine the gyroid 
structure, which was recently suggested as a different type of nuclear pasta by analogy with the 
polymer systems. The gyroid structure investigated in this paper is approximately formulated as 
' an extension of the periodic minimal surface whose mean curvature vanishes. In contrast to our 

^SJ , expectations, we find from the present approximate formulation that the curvature corrections act 

to slightly disfavor the appearance of the gyroid structure. By comparing the energy corrections 
^ ' in the gyroid phase and the hypothetical phases composed of d-dimensional spheres, where d is a 

general dimensionality, we show that the gyroid is unlikely to belong to a family of the generalized 
■ dimensional spheres. 

PACS numbers: 21.65.-f, 26.50.+X, 26.60.Gj 



I. INTRODUCTION 



Nuclei in laboratories are roughly spherical except for some significantly deformed nuclei, which are often encoun- 
tered in excited states. It is considered, however, that nuclei in matter just below normal nuclear density are not 
1 always spherical even in the ground state. Matter in supernova cores and in neutron star crusts are the best candidates 
\jQ ' where non-excited nuclei that have unusual shapes could be found, surrounded by a gas of dripped neutrons. The idea 
of non-spherical nuclei in matter was initially introduced when one considered a melting transition of matter with 
OO nuclei into uniform nuclear matter as the density increases. It was suggested that before melting, the system turns 
, into a bubble state in which nuclei turn inside out [H, Q • In the subsequent works using liquid drop models (LDMs) , it 
I ■ was concluded that nuclei deform from sphere (SP) to cylinder (C), slab (S), cylindrical hole (CH), and spherical hole 
' (SH) before melting into uniform matter [^-[sj. Since such non-spherical nuclei look like spaghetti, lasagna, macaroni, 
- and Swiss cheese, respectively, they are often referred to as nuclear pasta (for reviews, see, e.g., Refs. p-Q)- Recently, 
' by analogy with the polymer system, it was suggested that nuclei with a complicated shape called gyroid (G) may 
I> , appear between the C and S phases, and similarly, the hole structure of gyroid (GH) may appear between the S and 

CH phases Q . The nuclear shapes mentioned here are summarized in Fig. [T] 
rS , The possible presence of pasta nuclei in supernova cores and neutron star crusts would play an important role in some 
' astrophysical phenomena. For instance, pulsar glitches, which are thought to be caused by catastrophic unpinning 
of quantized neutron vortices in the inner crust of a neutron star, could be strongly affected by the configuration of 
nuclei in an amorphous solid that constitutes the crust [lo'| . It remains to be clarified how neutron vortices could pin 
in such an amorphous solid, particularly in the presence of pasta nuclei.^ In the context of core collapse supernova 
explosions and protoneutron star cooling, which are mainly controlled by the neutrino opacity, the possible influence of 
pasta nuclei on the neutrino opacity needs to be taken seriously. The cross section for the neutrino-nucleus scattering 
depends on the structure of nuclei d3- 15]. The neutrino scattering processes are no longer coherent in the directions 
in which nonspherical nuclei are elongated. It is different from the case of roughly spherical nuclei whose finiteness in 
any direction leads to constructive interference in the scattering. 

Earlier investigations of pasta nuclei basically focus on the zero-temperature phase diagram including the SP, C, S, 
CH, SH, and uniform phases. The Coulomb and surface energies are calculated in the LDM approach because their 
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^ Possible relevance of pasta nuclei to the vortex pinning was first proposed in Ref. in terms of pinning-induced nuclear rods. 
Subsequently, Jones |12|) concluded from the argument of formation enthalpy that pinning-induced nuclear rods are unlikely to form. 
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FIG. 1: Candidates for equilibrium nuclear shape. In this figure, the following notation is used: gyroid (G), gyroid hole (GH), 
sphere (SP), cylinder (C), slab (S), cylindrical hole (CH), and spherical hole (SH). 



delicate competition is responsible for the presence of pasta nuclei at sufficiently small internuclear spacings (e.g., 
Ref. [5]). In particular, the pasta phases inevitably appear between the SP and uniform phases at zero temperature 
in the incompressible limit where the density in the nucleus is fixed at the saturation density. Nevertheless, whether 
pasta nuclei actually occur in stellar matter remains to be answered. A compressible type of LDMs in which matter 
in the nucleus is allowed to have various densities helps to answer this question. Many calculations based on the 
compressible LDMs were performed by calculating the bulk and surface energies from a single model for the effective 
nucleon-nucleon interactions (e.g., Ref. [16]) and gave different predictions of the density region of pasta nuclei in 
neutron star crusts. In fact, such differences are shown to come mainly from the still uncertain density dependence 
of the symmetry energy [i7|. In the case of supernova cores where the system is less neutronrich but hotter than 
matter in neutron star crusts, the density region of pasta nuclei is predicted to be relatively wide at sufficiently low 
temperatures, while the description of the temperature-induced melting of pasta nuclei is still model dependent. 

In calculating the electrostatic energy of matter containing pasta nuclei within the LDMs, one typically assumes a 
geometrical shape of nuclei and uses the one-dimensional Wigner-Seitz (WS) approximation. For further calculations 
on the structure of pasta nuclei and a gas of dripped neutrons, one often goes beyond a simple liquid-drop picture 
and utilizes the Hartree-Fock and Thomas-Fermi methods within the WS approximation (e.g . , R efs. }lS-20j|). Beyond 
the WS approximation, shell effects on bound and unbound nucleons were considered [2l| - |23| : the shell effects on 
unbound nucleons are often referred to as band effects. On the other hand, multidimensional computations without 
any assumption of geometrical shape were also per formed within the framework of the Thomas- Fermi method [2^ , the 
quantum molecular dynamics (QMD) method |25l - l27| . and the Hartree-Fock theory (22i. [28l - [30j . which reproduce the 
five shapes listed above. Moreover, some of these multidimensional computations suggest more complex structures. 
We remark that pairing effects [l^ , electron-screening effects [HI, , and fluctuation-induced displacements of pasta 
nuclei [ssl - lssj have also been taken into account in describing the structure of pasta nuclei. 

It is noteworthy that similar phase diagrams are obtained for nanostructures of block copolymers (36| . Block 
copolymers are made up of two joined chemically distinct polymer blocks and each polymer block consists of a linear 
series of identical monomers. Since the same polymer blocks assemble spontaneously, phase separation of the domain 
structure occurs. The domains have a proper size and can not grow larger and larger because the distinct polymer 
blocks are joined strongly in each block copolymer. A resulting rubberlike entropic restoring force is counterbalanced 
by the thermodynamic interfacial tension driving macroscopic phase separation. This is a great contrast to the case of 
pasta nuclei where the Coulomb repulsion is essential to clustering with proper sizes as will be shown later. The shape 
of polymer domains depends on the fraction of two blocks while the shape of pasta nuclei depends on the volume 
fraction of the nuclear matter part. 

One of the shapes of polymer domains is the gyroid, which is a periodic bicontinuous morphology discovered 
experimentally for the regions between the C and S counterparts and the S and CH counterparts. Interestingly, in 
the case of nuclear pasta, some phases of complicated structure are predicted to appear in the corresponding regions 
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by the QMD simulations and Hartree-Fock calculations. Matsuzaki [37| pointed out the possibility that pasta nuclei 
may have such morphologies as observed in the polymer systems^ but made estimates that were too rough to address 
the tiny energy differences between morphologies. According to the precise Coulomb and surface energy evaluations 
within the LDM, it was found that the G (GH) phase does not appear in the ground state for any density but the 
energy difference from the most stable phase becomes quite small near the transition point from the C (S) phase to 
the S (CH) phase Furthermore, the volume fraction of nuclei at this point was shown to be ^ 0.35 (0.65), and 
this value is very close to the two-block fraction where the corresponding transition occurs in the polymer systems. 

The LDM approach to inhomogeneous matter at subnuclear densities was frequently utilized since the earliest studies 
■ It has advantages over other approaches not only because calculations of various thermodynamic quantities can 
be performed more straightforwardly, but also because various corrections can be added to a Weizsacker-Bethe mass 
formula in a systematic manner. One such correction is the curvature correction to the surface energy. The curvature 
effects, which are implicitly taken into account in the Thomas- Fermi and Hartree-Fock approaches, are often ignored 
in the LDM approach, but can be taken into account once the curvature coefficient is given. The importance of the 
curvature energy for nonspherical nuclei was pointed out by Pethick et al. 38]. They introduced this term so as to 
explain the difference of the nucleus-bubble transition density between Hartree-Fock and LDM calculations. They 
noted that, since the sign of the curvature corrections for hole nuclei becomes negative, the curvature corrections 
act to destabilize the phase with nuclei as compared with that with holes. In Ref. [39[, the curvature coefRcient was 
obtained by calculating the curvature thermodynamic potential in the Thomas-Fermi approximation with Skyrme 
interactions in a manner consistent with calculations of the surface tension. It was also shown that there are large 
variations in the curvature energy among the results for different interactions. Incidentally, the curvature term of the 
LDM is poorly known even from empirical nuclear masses I40l l4l| . The curvature corrections for other pasta phases 
were also investigated by utilizing the Hartree-Fock theory |42| and the Thomas- Fermi approximation 43] . 

In this paper, we build the curvature correction into the energy of pasta nuclei including the gyroid structure as 
derived within the LDM in Ref. |9|] and examine the associated change in the transition densities between different 
shapes. In this paper, we will also address the question of whether or not the gyroid phases are energetically favored 
by newly allowing for curvature effects. In addition, we will present details of the calculations of the surface and 
Coulomb terms omitted for want of space in Ref. Q . 

It is also interesting to examine the possibility that the transitions of pasta configurations are smooth as proposed by 
Ravenhall et al. f^. They regarded sphere, cylinder and slab as three-, two- and one-dimensional spheres, respectively, 
and extended the surface and Coulomb energy expression derived for these three specific cases to a geometry of general 
dimensionality d. By assuming cZ to be a continuous variable, they obtained a phase diagram in which the optimal 
value of d changes continuously with density up to a melting point into uniform matter. In this paper, we also discuss 
whether the gyroid, the recently proposed nuclear pasta, can be interpreted as a non-integer-dimensional sphere by 
comparing the Coulomb, the surface, and the curvature energies of the gyroid phase with those of a geometry of 
general dimensionality d. 

This paper is organized as follows. In Sec. [Hi we write down the expressions for the energy of pasta nuclei, which 
include the curvature correction. Section IIIII is devoted to illustration of the phase diagram of nuclear pasta, which 
allows for the possible gyroid phase, and to energy comparison between the gyroid and the general dimensional spheres. 
Our conclusions are presented in Sec. IIVI 

II. ENERGY OF MATTER AT SUBNUCLEAR DENSITIES 

In this section, we summarize the energy expressions that may be obtained for matter at subnuclear densities within 
the framework of the LDM, which include curvature corrections. We then derive the size equilibrium condition from 
the energy minimization by treating the curvature term perturbatively. 

A. Liquid drop model 

In describing the energy of matter at subnuclear densities, we regard nuclei as liquid drops containing neutrons 
and protons. Here, we deal with the zero-temperature matter, which consists of nuclei with a single shape and size, 



^ Instead of the gyroid morphology, he studied the double-diamond morphology, which is another periodic bicontinuous structure. His- 
torically, the double-diamond morphology was initially considered to be the most likely periodic bicontinuous structure for the polymer 
system. 
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FIG. 2: (Color online) A schematic density profile in a unit cell as described by the present LDM. The shaded areas in Fig. [T] 
correspond to "inside nucleus" of this figure. 



while nuclei with various shapes and sizes may be mixed at finite temperatures. In our model, nuclei occupy the 
volume fraction u and have the sharp boundary, uniform number density n'", and uniform proton fraction a;'". They 
are embedded in a uniform neutralizing background of electrons of number density Ma;'"n™ and, if any, in a gas of 
dripped neutrons of uniform number density ri°"*. For simplicity, we do not take surface diffuseness or neutron skin 
into account. Since the internuclear Coulomb repulsion gives rise to a spatially periodic structure, we consider a unit 
cell of volume by setting the shapes of nuclei a priori. In Fig. [2l a schematic density profile in a unit cell is shown. 

B. Energy expression including curvature term 

We write the energy expression in the form of the semiempirical mass formula which corresponds to the present 
liquid-drop picture. Including the curvature term and not including rest mass, the total energy of a nucleus with Z 
protons and A ~ Z neutrons can be written as 

-E'tot = + Eg + Ecuiv + Egym + -E-ssym + ^-Coul 

= (a„ + a,yn,(l - 2xf) A + {a, + a,,y„,{l ~ 2xf) A^/^ + a.^rvA^^^ + acoui^'^A^/^ , (1) 

where x = Z/A. The subscripts v, s, curv, sym, ssym, and Coul stand for the volume, the surface, the curvature, the 
symmetry, the surface-symmetry, and the Coulomb terms, respectively. This expression for the curvature term is also 
seen in Eq. (1.1) of Ref. [s^. Correspondingly, we write the total energy of a unit cell W as 

W^Wb+Ws+ Weurv + VKcoul , (2) 

where VFb, Ws, Wcurv, and IVcoui are the bulk energy, the surface energy, the curvature energy, and the Coulomb 
energy, respectively. This expression is the same as Eq. (3) of Ref. @ except that the curvature energy is introduced. 
Because of the saturation property of nuclear matter, the mass number A is roughly proportional to the volume of a 
nucleus. With the volume fraction of a nucleus, w, fixed, the volume of a nucleus is proportional to that of the unit 
cell a'^. Therefore, as we will see below, each energy term in Eq. ([2]) can be characterized by the a dependence. 
Since the bulk energy corresponds to the volume term in Eq. ([Ij, we can write it as 

= u;b(w,x'^n"^n°"')a^ (3) 

with the average bulk energy density x'", n™, n°"*). Note that we encapsulate all the energies that are propor- 

tional to the cell volume, such as the electron kinetic energy inside and outside the nucleus, in this term. 

The surface energy is proportional to the area of the surface of the nucleus, namely, the interface between a part 
of nuclear matter and a part of neutron matter. Thus it depends on the shape of the interface, whereas Wb does not. 
We write the surface energy of a cell as 

Ws = <j{x"\n"\ n°''')g{u, shape)a^ (4) 
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where 17(0;™, n™, is the surface tension, and the relative surface area g{u, shape) is the area of the surface for 
a = 1 with shape = SP, C, S, CH, SH, G, GH. Generally, it is given by 

g {u, shape) — — f dS, (5) 

J5(ti, shape) 

where the surface integral is taken over the nuclear surface S'(m, shape) within a unit cell. 
We turn to the curvature energy, which may be expressed as 

Wcurv = n'", shape)a, (6) 

where a;(a;'", n'", is the curvature coefficient. h{u, shape) is the integrated mean curvature defined as 



h{u, shape) = — / H{x,y, 

^ JS(«, shape) 



z)dS, (7) 



where H{x, y, z) is the mean curvature at the point (x, y, z) on S{u, shape). The definition of H is found in textbooks 
on differential geometry (e.g., Ref. [iJl) [see Eq. for the definition by implicit functions]. Note that shape) 
does not depend on a because H{x, y, z) is proportional to 1/a and is integrated over the nuclear surface. 
The Coulomb energy of a cell is written as 

WcouX = \ ! [e{np(r)-ne}0(r)]dr, (8) 

where e is the elementary charge. np{r) and rie denote the number densities of protons and electrons, respectively. 
In our model, their difference is expressed as 

, . 1(1 — u)a;'"n™, inside the nucleus, 

n„{r) — Uf, = < . . (9) 
I — ux'"/!'", outside the nucleus. 

4>{r) is the Coulomb potential determined by the Poisson equation: 

V20(r) = -47re [?ip(r) - rie] . (10) 



We can rewrite Eq. ^ as 



and Eq. (HU]) as 



VFcoui = (ea;'"n'") wcoui(w, shape)a^ (11) 



Vlqis) - -27rp{s). (12) 
Here, s — r/a and Vg are the dimcnsionlcss vector and its differential operator, 

wcoui(m, shape) = / p{s)q{s)ds, (13a) 

Jccll' 



p{s) 



[np{as) 



(13b) 



<?(«) = 7: — — T0(a*), (13c) 

where "cell"' in the integral denotes the cell normalized by a^. We can see that W^coui(u, shape) is proportional to 
a^, which is consistent with Eq. ((T|). 
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C. Energy minimization 

By substituting Eqs. ([3]), (HJ, ([6]), and (fTTj) into Eq. ([2]), we rewrite the total energy density as 

W , in in ont^ , (T(a:'" , n""^ )g(^, shape) ^ n'", shape) 

^ = Wb(u,x ,n ,n )H \ ^ 

a-^ a 

+ (ex'"n'")^ wcoui(u,shape)a^. (14) 

In the following, we minimize the total energy density. This consists of two steps. The first step is minimization with 
respect to the size of the unit cell a, which leads to 

Since the curvature energy is generally small compared with the surface and Coulomb energies, at first, we solve the 
equation for w = as 

/ \ 1/3 



For Luh <C <jgaQ, Eq. (flS)) can be solved as 



2(ea;'"n'")2?i;coui 



2ujh \ , , 

«o 1 + • (17) 



By eliminating a from Eq. (jl4p and by retaining terms of up to first order in tuh/agao, we obtain 
W 3 . i,, ^2/3 o/Q 1/3 / 2ujh \ 

73 - -^ + ^(-'"-'"-)' ^^'Voui(l + 5^j 

1/3 2/3 



^2/3 • 



(18) 



We remark that we can rewrite Eq. ([T5|) as Ws + 2Wcurv — 2T4^coui, an extended version of the well-known condition 
for size equilibrium, Ws — 2Wcoui, derived for w = 0. 

Second, we minimize Eq. (jl8p with respect to the shape for given u. This is performed by simply comparing the 
energies for different shapes and by finding the shape that gives the lowest energy. Note that, in Eq. (fT8|) . the shape 
dependence is entirely encapsulated in the geometrical factors defined as 

i^o(w, shape) — shape)^/"^u;coui(u, shape) (19a) 

u;coui(u,shape)2/^/i(M, shape) 

Fi u, shape = --j- . 19b) 

g(u, shape)"^/3 

Fq(u, shape) corresponds to the sum of the relative Coulomb and surface energy densities, which was denoted in our 
previous study [9] as -F(u, shape), whereas -Fi (u, shape) corresponds to the relative curvature correction. For a; = 0, 
as already discussed in Ref. [9|, the nuclear shape that minimizes the total energy density is determined uniquely for 
given volume fraction u and is independent of the average bulk energy density Wb and the surface tension a. Likewise, 
the shape dependence of the relative curvature correction is determined uniquely for given u. Because of the different 
shape dependences of Fq and -Fi, however, the nuclear shape is no longer determined uniquely for given u in the 
presence of the relative curvature correction. 

We remark that we will set specific values for the parameters n°"', n'", and a;™ by hand, although they can, in 
principle, be determined by additional energy minimization under boundary conditions associated with electric charge 
and baryon number. The setting will be made in the spirit of the incompressible LDM. 

III. EQUILIBRIUM NUCLEAR SHAPES 

In this section, we seek the equilibrium nuclear shapes by comparing the energy densities of the pasta phases, which 
include the gyroid phase. In particular, we address the question of whether the effect of the curvature is helpful for 
the gyroid appearance or not. We also compare the energy of the gyroid phase with that of a generalized dimensional 
sphere. 
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A. Candidate shapes 

As usual, we consider sphere, cylinder, slab, cylindrical hole, and spherical hole as candidates for the equilibrium 
nuclear shape. These shapes are conventional pastas and are regarded as the one-, two-, or three-dimensional spheres 
that minimize the surface area at constant nuclear volume. In addition, the gyroid phase, the nuclear pasta with a 
periodic bicontinuous morphology recently proposed by analogy with the nanostructures of block copolymers Q , and 
its hole structure are included in our investigations. 

For the conventional pasta phases, the relative Coulomb and surface energy densities can be calculated as 



Fo(u,SP) = (367 



2\2/9 



5?1 



-I 1/3 



+ CbccW 



(20a) 



(U - 1 - log U) + ChcxU^ 



1/3 



(20b) 



i^o(^,S) 



1/3 



,2/3 



(1 



N2/3 



(20c) 



The numerically determined coefficients, Cbcc = 6.5620 x 10^ and Chex — 1-2475 x 10^ , are corrections to the 
WS approximation, and the subscripts bcc and hex represent the body-centered cubic (bcc) and hexagonal (hex) 
lattices, respectively, which are the most stable alignments for each nuclear shape For an SHI and a CH, the 
following relations hold: Fo(u, CH) = Fq{\ — u, C) and Fo{u, SH) — Fo{l — u, SP) for the reasons as described below. 
First, the hole morphology with a fraction of u has the same surface area as the normal morphology with a fraction 
of 1 — w. Thus, the relations such as g{u, SH) = 5(1 — u, SP) hold. Second, since the charge density (|9]) has a 
sign opposite to the case of its hole morphology with u, the signs of p{s) defined by Eq. (|13bj) and q{s) defined by 
Eq. (|13cp also become opposite. According to Eq. (|13ap . therefore, the relative Coulomb energy satisfies the relation 

WCoul(l - SP) = WCoul{u, SH). 

Similarly, the relative curvature corrections can be obtained as 



Fi{u,SP) = 



V 3w 



1/9 



^ (2^5/3 _ + ^8/3) 
5^1 



ChccU 



2/3 



(21a) 



^^i(.,c)H^ 



1/3 



-1 2/3 



Y (U - 1 - log u) + ChexU^ 



(21b) 



Fi(u,S) = 0. 



(21c) 



For an SH and a CH, Fi {u, CH) = — Fi (1 — u, C) and Fi (u, SH) — —Fi{l — u, SP) are satisfied for the following reason. 
While the surface area of the hole morphology with a fraction of u is the same as that of the normal morphology with 
a fraction of 1 — u, the radii of curvature of the nuclear matter surfaces have signs opposite to each other. Thus, the 
relations such as h{u, SH) = —h{l — u, SP) hold. Incidentally, we remark that the volume fraction of the nucleus in a 
unit cell u has a geometrically allowed range: u < -s/Stt/S for an SP with bcc and u < tt/2\/3 for a C with hex. 

The gyroid structure is so complicated that we should evaluate g{u,G), /i(u, G), and uicouUw, G) numerically. We 
set the gyroid structure by the following level surfaces: 



, , . 2'KX 27rw 27rj/ 2t:z 27rz 2ttx 

J [X, y, z) — sm cos h sm cos h sm cos 

a a a a a a 



(22) 



where (x, y, z) are the spatial coordinates. Here, by taking a as the periodic length, the volume of a unit cube becomes 
a'^. By using this expression, we can assume that the region that satisfies \ f{x, y,z)\ > k corresponds to the nucleus in 
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FIG. 3: (Color online) (a) Relative Coulomb and surface energy densities, (b) their differences from the value of the slab phase, 
and (c) relative curvature corrections, calculated as functions of u for the seven phases of interest here. The notations are the 
same as in Fig. [T] 



Fig. [5J where fc is a positive parameter that specifies the volume fraction u of the nucleus in the unit cube.^ Note that 
fc — > corresponds to u — > 1 and that u is a monotonically decreasing function of k. One of the notable characters of 
this structure is that the regions inside and outside the nucleus are bicontinuous. Equation (|22p is no longer a good 
approximation for such small values of it as m < 0.0354 since the resultant configurations are not bicontinuous but 
pinched off. Because this poses no problem for the following analysis, we do not consider these configurations. In 
addition to the gyroid morphology, the hole structure of the gyroid, for which nucleons reside in the region satisfying 
\f{x, y, z)\ < fc, is taken into account. Again, the configurations are bicontinuous only for u < 0.965. 

To evaluate G) and h(u,G), we should perform the surface integrals ([5]) and ([7]) numerically. The mean 
curvature of a level surface defined by such an implicit function as Eq. (j22p is known to be given by 

rr/ \ I (.fxx ^" fyy)fz ^ (fyy ^ fzz)fx ^" ifzz ^ fxx^fy '^{fxyfxfy ^" fyzfyfz ^" fzxfzfx) 

^(->^.^) = ± ^UITf[TflW-^ ' ^''^ 

where f^y = gfg^ 2/? z) and so on (43 |. The sign ± corresponds to the normal and hole morphologies. Fortunately, 
in the cases of interest here, we can evaluate the surface integral by converting it to the volume integral via Gauss' 
divergence theorem. Mathematical details are given in the Appendix. This conversion is useful because the surface 
integral, in general, is difficult to perform numerically. 

We then evaluate wcoui(u, G) by solving the normalized Poisson equation ([T2|) . Here, we utilize the discrete Fourier 
transform, which is known to be powerful for a periodic cubic box. 



B. Curvature effect 



In Fig.|3l we show ^o(u, shape) and shape) as functions of m for all the shapes considered here. The differences 
from the value of the S phase are also shown for ^o(w, shape). For all the phases except the SP and C phases, Fq 
and Fi converge to zero at the point m = 1, which corresponds to uniform matter. Recall that for w = 0, the shape 
that minimizes F(){u, shape) is the equilibrium one. Thus, we find that equilibrium nuclei deform from SP to C, S, 
CH, and SH with increasing u. While the G phase does not give the minimum value of Fq for any m, the Fq value 
of the G phase is very close to those of the C and S phases at the transition point from the C phase to the S phase 



^ We remark that the surface with fc = is not a minimal surface itself, whose mean curvature vanishes but very close to it. This is 
because Eq. ((22} is just an approximation to the mathematical expression for the gyroid that is originally defined as a family of periodic 
minimal surfaces. For u < 1, the equilibrium configurations are most likely characterized by constant but nonzero mean curvature 
surfaces, which are a more general class of minimal surfaces that are stationary with respect to variations of the surface area for a fixed 
volume fraction. The mean curvatures of the surfaces l|22|l with corresponding k are not constant; they fluctuate by '^20% at m = 0.35. 
Again, this is caused by the approximation. Nevertheless, we call the surfaces I I22II "gyroid" in this paper. 
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{u = 0.35). The same is true of the Fq value of the GH phase at the transition point from the S phase to the CH 
phase {u — 0.65). 

The relative curvature correction is the largest for the SP phase, followed in order by the C, G, S, GH, CH, and 
SH phases. This hierarchy is due to the following reason. The curvature energy is a correction to the surface energy. 
Phcnomenologically, the surface energy arises because nucleons near the surface have less neighbors to interact with 
attractively than nucleons near the center. Therefore, the more curved surface inward, the less neighbors exist for 
nucleons near the surface. This feature indicates that the SP morphology has the largest mean curvature among the 
phases of interest here for the same volume fraction. The curvature of the S morphology is zero because of its flat 
surface. Since the surface of hole nuclei is reentrant, the sign of their curvature corrections becomes negative. 

To determine the nuclear shape that minimizes the total energy density (|18p at given u, we specify the values of the 
parameters n'", and x™, which, in turn, are related to the coefficients Wb, cr, and w. By bearing the application 
to matter in supernova cores in mind, we follow the setup in the incompressible limit as employed in Ref. 7] , while the 
formulations and analyses shown above are also applicable to the compressible case. We set — and n'" = uq, 
where ng = 0.165 fm~'^ is the saturation density. Then, the volume fraction is given by u — n/uQ, where n is the 
average nucleon number density. The proton fraction is set to x™ = 0.3. The surface tension is assumed to be 
a — 0.73 MeV fm~^, which comes from the liquid-drop parameters Ug = AnRQao and assym = — 47ri?Q(ToCsym with 
CTo = 1 MeV fm"^, Cgym = 1-7, and Rq = (3/47rno)^/'^, which reproduce the properties of isolated finite nuclei in the 
limit of M — ?> 0. We assume that the curvature coefficient cj is a constant free parameter because of its uncertainties. 
In Fig. 21 we show the average total energy densities minus the bulk term for several values of a;, which are taken 
in the range w < 1 MeV fm~^ by reference to the calculations from the Skyrme interactions [3^ |4^. We can 
recognize that the curvature effect on this energy density difference is at most ~10% and that the assumption that 
the curvature energy is small is good for a; < 1 MeV fm~^. Note that we do not have to take a specific value of wi, in 
the incompressible case as considered here, but Wb is essential for the realistic description of the melting into uniform 
matter [l7j . 

In Fig. [SJ the difference AW/a^ between the average total energy density for each of the shapes considered here 
and that of the most stable phase at given n/ng is shown for various values of w. Whereas we originally expected 
the G phase to appear as the ground state because of the curvature correction, this is not the case. In the absence 
of the curvature correction, the total average energy density of the G morphology becomes the same as that of the 
C morphology at n/no ~ 0.35. The density where the C and G morphologies have the same energy is lowered by 
the curvature correction because the mean curvature of the G morphology is smaller than that of the C morphology. 
Therefore, the G phase is expected to appear at n/no ^ 0.35. However, the energy of the S morphology becomes lower 
than that of the G morphology because the curvature correction acts to raise the energy of the G morphology while 
keeping that of the S morphology unchanged, which has zero curvature. Thus, the G morphology does not appear 
for any density. In Fig. [51 we show the phase diagram on the n/no versus the w plane. We can recognize that, for a 
larger curvature coefficient, the transition densities between different morphologies become lower. Again, the reason 
is that the mean curvature of the low-density phase is larger than that of the high-density phase. 
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0.2 0.4 0.6 0.8 1 
n / 

FIG. 6: Phase diagram on the n/no versus the lj plane, calculated from the LDM with the parameters used in Ref. [3]. The 
notations are the same as in Fig. [1] 



C. General dimensionality 

In the earliest study on the nuclear pasta [1] , the Coulomb plus surface energy of the d-dimensional sphere was 
represented by a single expression, and it was argued that the continuous change of the nuclear shape with density 
could be characterized by the optimal value of the continuous variable d. Here, we investigate whether the gyroid can 
be regarded as a member of the d-dimensional spheres or not. In this subsection, we neglect the corrections to the 
WS approximation for simplicity. By following a line of argument from Ref. , we write the relative Coulomb and 
surface energy density Fq and the relative curvature correction Fi for the d-dimensional spheres as 



Fo{u,d) 



d + 2 



du^-y^' - 
d-2 



1/3 



(24a) 
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FIG. 7: (Color online) Same as Fig. [3] but for the d-dimensional spheres and gyroid. 



d-1 



2du 



2{d + 2y 



d-2 



21 1/3 



Note that, 



lim 



du^-^/'' - 
d-2 



= logu + 1, 



(24b) 



(25) 



and d = 3, 2, and 1 correspond to sphere, cylinder, and slab, respectively. 

In Fig. [71 we compare Fq and Fi of the gyroid with those of the c?-dimensional spheres with 1 < d < 2. Least 
squares fitting allows us to determine the "dimension" of the gyroid as c? = 1.478 for Fq and d = 1.515 for Fi. In 
this sense, the G phase could be interpreted as the intermediate of the C and S phases. As can be seen in the figure, 
however, the general dimensionality does not give a good overall fit for Fq and Fi of the gyroid. This fact suggests 
that the candidates for equilibrium nuclear shape should be examined individually beyond the general dimensionality. 



IV. CONCLUSIONS 



In this paper, we built tlie curvature correction into the LDM and examined its influence on the equilibrium phase 
diagram associated with the nuclear pasta including the gyroid structure. We confirmed that the curvature effect 
pushes the onsets of the shape changes to low densities. We originally expected the gyroid to appear in the phase 
diagram because of the curvature correction. However, we found that the curvature correction makes the gyroid 
appearance harder although the effect is not remarkable. We also compared the energy of the gyroid with that of a 
generalized dimensional sphere. We found that the G phase does not belong to a family of such spheres. This fact 
implies that the intermediate phases between the conventional pastas should be described beyond the scope of the 
noninteger dimensionality to know their possible existence in the ground state. 

While we mainly considered supernova matter, the qualitative results for the curvature effect would be unchanged 
for neutron star matter. In this case, the proton fraction inside a nucleus, a;'", is smaller in the absence of trapped 
neutrinos, and a gas of dripped neutrons appear in the inner crust of a neutron star. Accordingly, the coefficients 
Wb, (T, and uj are greatly modified while Fo(u, shape) and Fi(m, shape), Eq. (|19l) . are unchanged. In the present 
LDM approach, we concluded that the G phase docs not appear in the equilibrium phase diagram. However, the 
energy difference between the G phase and the ground state is so tiny that it would be significant to confirm this 
conclusion in more sophisticated approaches. We itemize various corrections ignored in the present analysis as follows: 
(i) Compressibility. It leads to deviation of the liquid drop density from uq even in equilibrium, which in turn results 
in changes in the bulk energy (|3]), the surface energy (|3]), the curvature energy ([5]), and the Coulomb energy ([S]). To 
deal with the compressible case, we should assume a specific nucleon-nucleon effective interaction, which gives the 
density dependence of Eqs. (01), (H]), and (|5]). Also, note that compressibility leads to surface diffuseness and 
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charge screening as categorized below, (ii) Surface diffuseness. The nuclear surface is not a sharp boundary because 
of finite compressibility. The resultant surface diffuseness corrects the Coulomb energy [l| . (iii) Neutron skin. The 
mean-square radius of the neutron distribution is generally larger than that of the proton distribution, which corrects 
the bulk energy ^. (iv) Charge screening. Both protons and electrons redistribute and then contribute to a 
decrease in the Coulomb energy (|8]) and an increase in the bulk energy ([3]) [l|, HH, [H, H^] . (v) Coulomb exchange 
energy. The proton Fock term corrects the Coulomb energy ( 8) 111 , (vi) Thermal fluctuations. Thermally induced 
displacements tend to destroy the nuclear pasta structure f33|-[35|. It is also interesting to ask whether or not the 
coexistence of the G phase with other phases is possible at finite temperatures. Note that the G phase tends to appear 
at finite temperatures for the copolymer systems, because the entropic restoring force, which drives microscopic phase 
separations, can be dominant over the destructive effect by thermally induced displacements. 

In the course of this work, we have devised a systematic method of calculating the Coulomb, surface, and curvature 
energies of nuclei of various shapes, which include the one that has periodic bicontinuous structure. In light of the fact 
that pasta nuclei of complicated shapes have recently been taken note of [1, [2^, [13] , we expect that this method 
could be useful for future theoretical work involved on pasta nuclei. 
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Appendix A: Calculation of the surface integrals 

Equations (O and ([7]) include the integrals over the nuclear surface in a unit cell, which are difficult to calculate 
directly for such complicated surfaces as the gyroid. In this Appendix, we show the numerical method to evaluate 
these integrals for the gyroid. As mentioned above, a unit cell of the gyroid is a periodic cubic box with volume a^. 
Thus, we consider a more general surface S'(u, shape) and an arbitrary continuously differentiable vector field v(r) 
that are defined in a unit cube and satisfy the periodic boundary condition. From Gauss' divergence theorem, one 
can write 

(j) v{r) ■ n{r)dS = / V • f (r)dr, (Al) 

<y S <y V (u,sha.pc) 

where n(r) is the outward pointing unit normal vector of the surface S and the volume integral is taken over the 
region, ^(m, shape), surrounded by the surface S. Note that the surface integral in Eq. (jAip is taken over the cell- 
nucleus boundaries plus S'(m, shape). For simplicity, we consider a case in which the integration region is periodically 
connected only at the top and bottom sides of the cell as shown in Fig. [51 Then the integral can be divided as 



* v{r)-n{r)dS= / v{r) ■ n{r)dS + / v{r) ■ n{r)dS + / v{r) ■ n{r)dS. 

Js JS(m, shape) Ja Jb 



(A2) 



The integrals over the top side A and bottom side B are canceled out because t;(r)'s on A and B are parallel because 
of the periodic condition but n(r)'s on A and B are obviously antiparallel. Consequently, only the integral over the 
nuclear surface shape) remains. Thus, we get 

/ v{r) ■ n{r)dS = / V • v{r)dr. (A3) 

^(ti, shape) shape) 

Substitution of v{r) ~ n{r) and H{r)n{r) into Eq. (IA3|) enables us to calculate the integrals in Eqs. ([5]) and (O, 
respectively. 

The unit normal vector of a level surface defined by such an implicit function as Eq. is ±V/(r)/|V/(r)|, where 
the sign ± corresponds to the inward and outward pointing vectors. Note that there are a finite number of singular 
points where V/(r) = 0. Fortunately, in the present case, these singular points are removable because the singular 
points are isolated ones rather than lines and planes. First, let us remove small spheres which enclose the singular 
points from the integration region of the right-hand side of Eq. (|A3[) . Correspondingly, the integrals over the surfaces 
of the small voids are added to the left-hand side of Eq. (jA3[) . These surface integrals vanish in the limit e — )■ 0, 
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S {u, shape) 




FIG. 8: (Color online) A schematic of a periodic nuclear surface in a unit cube. 



where e is the radius of the removed spheres. This is obvious in the case of v{r) = n{r) because the absolute values 
of the surface integrals are, at most, Aire^. In the case of v{r) — H{r)n{r), the surface integrals are 0{e). Thus, we 
can evaluate the surface integrals in Eqs. ([5|) and ([7]) by numerically computing the volume integrals that have the 
singular points removed. 
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